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1. INTRODUCTION 

Shock-boundary layer interaction can significantly influence not only the local 
transonic flow on missiles, wings and turbine blades but its influence can also ex¬ 
tend downstream within the boundary layer and thereby alter the global aerodynamic 
properties of lift, drag and pitching moment. It is therefore important that shock¬ 
boundary layer Interactions and their Reynolds and Mach number-scaling be properly 
modeled in engineering flow field prediction methods for supercritical aerodynamic 
bodies. This paper describes the application of a non-asymptotlc triple deck theory 
of transonic shock-turbulent boundary layer interaction which provides such a tool 
for non-separating two-dimensional flows over a wide range of practical Reynolds 
numbers. Section 2 contains a brief description of the essential features of the 
theoretical model. Sections 3 and 4 then describe how the results of a comprehen¬ 
sive parametric study of this theory may be used to develop a generalized "viscous 
wedge" model of the local interaction which embodies proper scaling behavior as 
well as an approximate account of incipient separation that is in good agreement 
with experimental trends. In Section 5 we examine the application of this theory 
as an element in global viscous flow field analyses of supercritical airfoils. In 
such problems it will be shown that even in non-separating cases the changes across 
the interaction may significantly alter the subsequent turbulent boundary layer be¬ 
havior for appreciable distances, especially when larger downstream adverse prersure 
gradients are present. 

2. BRIEF OUTLINE OF THE LOCAL INTERACTION THEORY 

Unlike significantly-separated flow where the disturbance flow pattern as-orI 
with a nearly-normal shock-boundary layer interaction is a very complicated an^ in¬ 
volving a bifurcated shock pattern', the unseparated case pertaining to turbulent 
boundary layers up to roughly Mj « 1.3 has instead a much simpler type of lp»era<. tiun 
pattern which is more amenable to analytical treatment. With some judic<ous simpli¬ 
fications, it is possible to construct a fundamentally-based approximate theory . f 
the problem in the latter case. Consider a known adiabatic boundary layet pro 11 ><s 
Mgfy) subjected to small transonic disturbances due to an impinging weak and nearly 
normal shock. In the practical Reynolds number range of interest here (10 5 <Re^< r), 
it has been established 2-4 that the local interaction disturbance field in the ne gh - 
borhood of the impinging shock organizes itself into three basic layered regions >-r 
"decks" (Fig. 1): 1) an outer region of potential invlscld flow above the boundary 
layer, rrtiich contains the incident shock and interactive wave systems; 2) an late' 
mediate deck of frozen shear stress-rotational lnviscid disturbance flow occupy;r.* 
the outer 90X or more of the incoming boundary layer thickness; 3) an inner shear 
disturbance sublayer adjacent to the wall which accounts for the interactive skit 
friction perturbations (and hence any possible incipient separation) plus most of 
the upstream Influence of the interaction. The "forcing function" of the problem 
here is thus impressed by the outer deck upon the boundary layer; the middle 
couples this to the response of the inner deck but in so doing can Itself modify 
the disturbance field to some extent, while the slow viscous flow in the thin Inrn 
deck reacts very strongly to the pressure gradient disturbances Imposed by these 
overlying decks. In treating this interactive field we employ a non-asymptotlc an*’.fod 
that is an extension to turbulent flow of Llghthill's approach 6 because of Its es¬ 
sential soundness and adaptability to practical engineering problems, similarity tc 
related types of multiple-deck approaches that have proven highly successful In treat¬ 
ing turbulent boundary layer response to strong adverse pressure gradients, and the 
large body of turbulent boundary layer Interaction data plus recent numerical studic- 
wlth the full Navler-Stokes equations which support the predicted results in a varte*v 
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of problems (see the survey in Ref. 5). Moreover, this approach provides at realistic 
Reynolds numbers a treatment of the inner deck pressure gradient terms plus the 
middle deck 3p/3y and streamline divergence effects, along vith simplifying ap¬ 
proximations that render the resulting theory tractable from an engineering stand¬ 
point . 

A very detailed description of the above-mentioned non-asymptotlc triple deck 
analysis can be found in Refs. 5 and 7 and hence will not be given here. The re¬ 
sulting predictions, such as typically illustrated in Fig. 2, describe all the es¬ 
sential global features of the mixed transonic character of the non-separating norma) 
shock-turbulent boundary layer interaction problem including the interactive pressure 
distribution and upstream influence, displacement thickness and local shape factor, 
and interactive skin friction up to incipient separation. This interaction theory 
employs for the Incoming turbulent boundary layer velocity profile a very general 
Composite Law of the Hall-Law of the Wake profile model due to Walz 8 , which is 
characterized by three parameters (Mi, boundary layer thickness Reynolds number and 
the incoming shape factor). The influence of both shock obliquity^ and wall curva¬ 
ture 10 have also been examined in detail and incorporated into the theory. Very 
extensive parametric studies and detailed comparisons with experiment have shown 
that it gives a very good account of all the important engineering features of the 
interaction over a wide range of Mach and Reynolds number conditions. Moreover, the 
Important but heretofore-ignored influence of incoming boundary layer shape factor 
Hj^ (and hence the upstream pressure gradient history) has been determined^* 7 . 

3. A GENERALIZED VISCOUS-RAMP MODEL OF THE INTERACTION 

In certain engineering applications to global flow field analysis computer o»> 
grams for wings or turbine blades, it has proved expedient to model the interaction 
as a simple local lnvlscld 6* - "bump" or "ramp." A serious deficiency of this ap¬ 
proach is that it does not account for the dependence of the bump shape aed size on 
Reynolds number, shock strength and boundary layer shape factor, while the additional 
interaction effects on the downstream boundary layer (such as Cf reduction; ire ig¬ 
nored altogether. With the aforementioned parametric study results in hand, hovpver, 
the present theory provides a much improved'Viscous ramp" representation of che in 
teraction (see Fig. 3A) whose key physical features have the correct dependence ci 
Mi, Re^ and 

Results for these viscous wedge properties taken from Ref. 11 are presented ;n 
Figs. 3B-3E, where the upstream and downstream influence distances, the slope and 
overall height of the 6* - bump are plotted along with the downstream C t /Cf Q values 
that may be needed along with 6* to re-initialize a subsequent turbulent boundary 
layer calculation downstream. We note in general that the overall scale of the to 
teraction (which can be a sensitive effect in both steady and unsteady applications 
where such viscous wedge models are employed) does not scale according to the un¬ 
disturbed boundary layer thickness even in the non-separating case. It is further 
noted that the viscous wedge slopes are in rough agreement with the maximum attached 
shock deflection value observed empirically 12 , although here of course there Is a 
dependence on Re^ and Hj. as well as Mach number. Finally, in all of these curves 
we see a significant dependence on the incoming boundary layer shape factor that wouJd 
appear to be an Important consideration in practical applications. 

Closed-form analytical fits to the various curves of Fig. 3 have been deveJope<> ‘' 
which provide a very rapid yet complete modeling of the local Interaction effects to- 
incorporation as a local module in global lnviscid-boundary layer analysis or design 
programs for supercritical airfoils. 
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4. INCIPIENT SEPARATION 

The present theory, although it breaks down at separation, does yield a useful 
indication of incipient separation where Cf . 0, owing to the particular attention 

paid to the treatment of the local interactive skin friction behavior • Since this 
indication is of great practical interest, a parametric study of incipient separa¬ 
tion conditions inherent in the present theory was carried out; the results for a 
normal shock on a flat surface are presented in Fig. 4A where the shock Mach number 
above which incipient separation occurs is plotted as a function of the Reynolds num¬ 
ber with the shape factor as a parameter; also shown in the Figure is the approximate 
experimental boundary determined by a careful examination 11 of a large number of trans¬ 
onic Interaction tests, besides Nussdorfer* s 1 ® M ~1.30 criterion for turbulent flow. It 
is seen that the theoretical prediction of a gradual increase in the incipient separa¬ 
tion Mach number value with Reynolds number is in agreement with the trend of the 
data; moreover, the theoretical prediction of only a small influence of shape factor 
on the incipient separation conditions is also borne out by the lack of any consistent 
H-effeet for the name Re discernible in the data (Squire has observed a similar insen¬ 
sitivity to Hj. in purely supersonic flow interactions 14 ). The absolute values of 
the incipient separation Mach number predicted by the present interaction theory are 
seen to be consistently slightly lower than the average experimental value; this is 
attributable to the combined effects of the linearized inner deck theory (which over- 
predicts the pressure gradient effect on Cf) and the assumption of a normal shock 
when in fact most of the experiments likely entail some shock obliquity (which also 
delays separation to somewhat higher shock-strengths). It is interesting to note 
that Nussdorfer's 13 original incipient separation criterion, based as it was on a 
very limited base, does roughly go through the average of the data although it does 
not account for the proper Reynolds number effect. 

Fig. 4B shows the influence of wail curvature; it is seen to have onlv a small 
effect in delaying incipient separation to a slightly higher shock strengtr> for » 
given Re^ and Hf, this being of the same order as the experimental data b„r.d. 

5. APPLICATION TO GLOBAL FLOW FIELD ANALYSIS 
5.1 Downstream Effects from a Shock-Boundary Layer Interacti on 

In addition to the Increased displacement thickness, the foregoing discussion 
shows that the skin friction level following the interaction is significantly redu.-d; 
combined with the attendant sensitivity to the profile shape, this suggests that 
subsequent downstream boundary layer development may retain a memory of the Intetac¬ 
tion effects for a considerable distance (over and above a simple thickening), par 
ticularly as regards possible incipient separation in any adverse pressure giadleo-s 
downstream. As indicated in Fig. 5 this downstream "interaction after-effect" in 
the boundary layer Influences the sensitive trailing edge region and thus may be ia- 
portent in the design and analysis of rear-loaded airfoils, especially at higher ltn 
coefficients with increasingly-aft shock locations; it likewise may be Important on 
three dimensional wing configurations where the shock interaction zones are well «fr. 

The aforementioned after-effect question was therefore subjected to detailed 
study using the two-layer turbulent boundary layer program of Moses ls as a model ct 
the downstream viscous flow; the program is coupled to the present interaction theot> 
by initializing it with the post-interactive flow properties so as to account fully 
(both Cf and 4*), partially (&* only) or not at all for the changes across the inter¬ 
action. Calculations were then made of the subsequent downstream turbulent bounder* 

behavior (8, C-, 8*, 6*) in various constant post-shock adverse pressure gradient* 
typical of airfoils for different assumed local interactive shock strengths and post 
tions or Reynolds numbers. The results serve as a paradigm of the downstream sensi¬ 
tivity question in real flows. 

A variety of cases were studied 16 , typical results of which are presented in 





Fig. 6 where we show Che predicted behavior of the boundary layer shape factor and 
skin friction in three lncreasingly-strong adverse pressure gradients downstream 
of an interaction occurring at a typically rearward position; the consequences of 
fully, partially or negligently-treating the boundary layer changes across the inter¬ 
action are indicated. Generally, it is seen that the downstream behavior of the 
boundary layer is indeed sensitive to detailed modeling of the interactive effects 
and that this sensitivity increases with the strength of the downstream adverse 
pressure gradient. The adverse pressure gradient magnifies the subsequent influ¬ 
ence of the skin friction (as well as the £*-rise) across the Interaction so that 
downstream separation tends to occur earlier than would be predicted by either neg¬ 
lecting or treating only the 6* effect of the upstream interaction. As shown in 
Fig. 7 , these predictions are supported by a comparison with boundary layer measure¬ 
ments downstream of a non-separating shock interaction zone on a supercritical airfoil; 
both the skin friction and shape factor data are poorly predicted when the interaction 
is neglected but are reasonably well predicted when the complete interaction effects 
are taken into account. 

Examination of many such results leads to the further conclusion that such in¬ 
teractive after-effects extend at least 20-302 chord distances downstream on a typical 
airfoil or wing and increase (as expected) with either larger shock strength or 
decreasing Reynolds number. If the trailing edge region lies within this range of 
the shock, it is thus seen that a simple thickening effect alone is not sufficient 
to account for the interaction and may result in an inaccurate prediction of the 
rearward boundary layer shape factor, skin friction and incipient separation proper¬ 
ties including their scaling. This is of practical importance for two major reasons: 

1) in regions of sustained adverse pressure gradient that often follow the short- 
scale interaction zone, the shape of the velocity profile and streamwise shear stress 
distribution (as well as thickness) are of considerable importance to the aerodynamic 
design of an airfoil or wing; 2) the altered boundary layer properties (especially 
possible incipient separation) near the trailing edge and into the wake can further 
exert a powerful effect on the overall aerodynamics via their influence of the Kutta 
condition 17 and on possible buffet onset. 

5.2 Supercritical Wing Section Flow Fields 

Nandanan et al 12 have carried out an even more detailed study of interactions on 
actual supercritical airfoils including experimental comparisons. They developed a 
global computational method for transonic airfoil flow analysis which incorporates the 
present analytical solution for near-normal shock-boundary layer Interaction into a 
state-of-the-art viscous-lnviscid computation code. Theoretical results obtalneo 
with this method were compared to representative data from boundary layer and surface 
pressure measurements on three transonic airfoils in the DFVLR-AVA (Gottingen) Trans¬ 
onic Hind Tunnel; some examples of these comparisons are shown in Figs. 8A and 8B The 
agreement between theory and experiment in both the boundary layer displacement chick 
ness and the surface pressure distributions was, for all test cases considered, quite 
good. The associated predictions of the local skin friction variation through the 
interaction zone also agree reasonably well with the values Inferred from the experi¬ 
mental boundary layer profiles via the Ludwig-Tlllman relation (see, e.g.. Fig. 9 ). 

The results of this investigation indicated that treating the shock-boundary layer 
interaction by conventional boundary layer theory generally leads to a slight underpve 
diction of the displacement thickness immediately downstream of the shock and, due ra 
the amplifying effect of the sustained rear adverse pressure gradients, to an appre¬ 
ciable underestimation of the displacement thickness at the trailing edge (see Fig.10). 
The latter is also clearly reflected in the pressure distributions and aerodynamic 
coefficients compared. Considering these results, one may conclude that it is general y 
necessary to include a physically correct treatment of shock wave-boundary layer intei 
action in the analysis of transonic airfoil flow. 
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6. CONCLUDING REMARKS 




The results of this study have shown that It Is now possible to Incorporate 
as an Interactive module within a global flow field analysis the correctly-modeled 
(and scaled) local shock-boundary layer interaction effects for the non-separating 
case. The non-asymptotic triple deck interaction theory involved covers a wide 
range of practical Reynolds numbers and turbulent boundary layer profile shape fac¬ 
tors including the effect of wall curvature; moreover, it gives an approximate in¬ 
dication of when incipient separation occurs. It was further shown that such theory 
is generally desirable when accurate predictions are desired in the important trail¬ 
ing edge region of rear-loaded supercritical airfoils because the detailed changes 
across an upstream interaction can significantly alter the subsequent turbulent 
boundary layer behavior for appreciable distances downstream. 
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